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id ■ ABSTRACT 

Cn I We present results from a monitoring program of 42 quasars from the Palomar- 

OV , Green sample. The objects were observed for seven years at the Wise Obser- 

0\ ■ 

■~-^ i vatory, as part of a long term effort to monitor AGN of various types. This is 

^ ■ 

f^ i' the most extensive program of its kind carried out to date on a well-defined 

O ' optically-selected quasar sample. The typical sampling interval is '^ 40 days. 

C/3 , One third of the quasars were observed at ~ 60 epochs and the rest at ~ 30 

• • I epochs in two bands {B and R) with photometric accuracy of ~ 0.01 mag. 
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We present lightcurves for all of the sources and discuss the sample variabil- 



Vh ' ity properties. All of the quasars in the sample varied during the campaign 

with intrinsic rms amplitudes of 5% < as < 34% and 4% < an < 26%. The 
rms amplitude and colour for the entire sample are ub = 14%, ur — 12%, 
and aB-R — 5%. On time scales of 100 to 1000 days the power spectra of 
the sources have a power-law shape, P^, oc v^'^ , with 7 « 2.0 and a spread 
^ 0.6. At least half of the quasars, particularly those that are most variable, 
become bluer when they brighten, and the rest do not show this behaviour. 
We quantify this phenomenon, which has been observed previously mainly 
in Seyfert galaxies. The quasars which are most variable tend also to exhibit 
asymmetry in their variations, in the sense that the brightening phases last 
longer than the fading phases. We have searched for correlations between the 
measured variability properties and other parameters of the quasars, such as 
luminosity, redshift, radio loudness, and X-ray slope. We find several new cor- 
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relations, and reproduce some of the correlations reported by previous studies. 
Among them are an anticorrelation of variability amplitude with luminosity, 
a trend of the autocorrelation time scale with luminosity, and an increase in 
variability amplitude with H/3 equivalent width. However, all of these trends 
have a large scatter despite the low observational uncertainties. 

Key words: galaxies: active - quasars: general. 



1 INTRODUCTION 

Continuum variability is a common property of active galactic nuclei (AGN), and is observed 
in all accessible spectral bands, from gamma-rays to radio. AGN liglitcurves show structure 
over a wide range of time scales and amplitudes. It is now well established that there are at 
least two variability classes of AGN. The first includes BL Lac objects and optically violent 
variable (OVV) quasars, whose variability is characterized by short time scales and large 
amplitudes - properties which are usually attributed to relativistic beaming (e.g., Bregman 
1990). The second class, which includes the majority of AGN (~ 90%), is the smaller- 
amplitude variables, which show variations of up to tens of percents in amplitude on time 
scales of weeks to years. The physical origin of their variability may be completely different, 
e.g., instabilities in accretion disks around massive black holes (Siemiginowska & Elvis 1997; 
Kawaguchi et al. 1998). 

Among the many studies of AGN variability in the various spectral bands, many workers 
have attempted to quantify the optical variability of quasars, and have discussed its correla- 
tions with redshift, luminosity, and spectral properties. Table |l| summarizes the observations 
and the results of some major studies of quasar variability. We list only the most recent 
paper from each program. 

The motivation for this extensive research has been the hope to constrain the nature of 
the central engines of quasars using variability data. Unfortunately, even basic results, such as 
whether luminosity and variability amplitude of quasars are correlated, are discordant among 
different studies. This suggests that sample selection, temporal sampling, and measurement 
errors may be affecting the conclusions. In early studies, samples were often chosen arbitrarily 
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Table 1: Summary of previous studies of quasar variability. Variability amplitude (cr) and variability timescale (x) ^vere defined in different 
manners by the various authors. Positive correlations are marked ^vith " -(-" and anticorrelations ^vith " — ". "0" means that the trend ^vas 
studied but no significant correlation ^vas found. 
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and concentrated on bright or "interesting" objects, that showed large variations. Such 
samples may be unrepresentative of the population of quasars as a whole, and the results 
may simply reflect the selection criteria. Even well-defined samples may produce spurious 
results due to irregular sampling. The sampling of most previous monitoring was sparse, 
making the variability properties sensitive to a few data points. The lightcurves were usually 
also unevenly sampled, complicating most time series analyses. 

We are engaged in several long-term quasar monitoring programs at Wise Observatory 
that attempt to overcome some of these problems. Netzer et al. (1996) described a 3-6 
year monitoring program of radio-loud quasars. The present paper presents the results of 
a program aimed at an optically selected quasar sample. We have monitored a well-defined 
subset of the Palomar-Green QSO sample (PG; Schmidt and Green 1983) for seven years. 
The PG sample is one of the best-studied optically-selected quasar samples. Its properties 
have been quantified at many wavelengths: radio - Kellermann et al. (1989), IR - Neugebauer 
et al. (1987), optical- Miller et al. (1992), and X-ray- Laoret al. (1997). The optical emission- 
lines and their variability in the PG sample were studied by Boroson & Green (1992), Maoz 
et al. (1994), and Kaspi et al. (1996), and the sample's X-ray variability properties were 
studied by Fiore et al. (1998). The polarization properties of the PG sample were studied 
by Berriman et al. (1990). 

Comparison between the various studies listed in Table |l] needs to be done with care, 
because the programs sample different regions of the absolute magnitude-redshift plane. In 
the present work, most of the sample are low-luminosity (and low-redshift) quasars (see 



§ |2.1| ). By the commonly-adopted absolute-magnitude criterion for quasars, Mb < —22, a 
minority of the sample may be considered Seyfert galaxies. Our sample and the samples 
of Bonoli et al. (1979), Giallongo et al. (1991), Cimatti et al. (1993), and Trevese et al. 
(1994) consist of quasars of similar absolute magnitudes (—28 < Mb < —21, see Figure P, 
but the redshift ranges in those studies are broader. Variability studies of bona-fide Seyfert 
galaxies have mainly focused on individual objects (e.g., Clavel et al. 1991), chosen by virtue 
of their high variability. The work of Bershady et al. (1998) is an attempt to quantify the 
variability of Seyfert galaxies in general. In terms of observational parameters, the sampling 
rate of our program is much higher and more even than that of past programs. The high 
accuracy attainable with differential CCD photometry allows us to probe smaller amplitude 
variations. 



© 1998 RAS, MNRAS 000, 



Long-term optical variability properties of the Palomar-Green quasars 5 

In §0 we describe the sample, the observations, and the reduction. In §|^ we describe the 
time-series and correlation analyses. In §^ we summarize our conclusions. 

2 OBSERVATIONS 

2.1 The Sample 

Our goal is to investigate the long-term variability properties of a particular sample of 
optically selected quasars and their connection to other quasar properties. The PG sample is 
well suited for this purpose, being a statistically complete sample (however, see Goldschmidt 
et al. 1992; Kohler et al. 1997) of UV-excess quasars. We chose from the PG sample all 
nearby (redshift z < 0.4), bright [B < 16 mag according to Schmidt & Green 1983) quasars 
(absolute B magnitude Mb < —22 mag, assuming a Hubble constant Hq = 70 km s~^Mpc~\ 
a deceleration parameter go = 0.2, and a power-law-shaped continuum, F^, oc z/~", with 
spectral index of a = 0.5 for the K-correction; see Schmidt & Green 1983) with northern 
declinations. We also required that there be a star within 4 arcminutes of the quasar with 
comparable brightness. The redshift criterion was applied so that the main optical emission 
lines (Ha and H/5) remain within the optical range for the purposes of a spectroscopic 
study (Maoz et al. 1994; Kaspi et al. 1996) carried out in parallel to the present work. The 
properties of the 42 PG quasars passing these criteria are summarized in Table 2 in the 
following format: 

Column (1) — Right ascension (1950). 

Column (2) — Declination (1950). 

Column (3) — Redshift. 

Column (4) — Number of observing epochs. 

Column (5) — Total observing period, in days. 

Column (6) — Median sampling interval, in days. 

Columns (7)- (10) — Median apparent and absolute magnitudes in the B and the R 

bands based on our measurements. The absolute magnitudes were computed assuming a 

spectral index, 

mediB~R)-{Bo-Ro) , ^ ..,, 

2.5 log ^ 

where Bq and Rq are the zero-points of the Johnson-Cousins magnitude system (Allen 1983), 

and Xb = 4400 A and Xr = 6400 A are the central wavelengths of the B and the R bands, 
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Figure 1. Distributions of absolute B magnitudes of the sample. The original distribution used for sample selection is repre- 
sented by empty bars, and the distribution based on our measurements is represented by filled bars. 

respectively. Using our measured colours and time-averaged magnitudes, the sample distri- 
bution of absolute magnitudes is somewhat changed from that used to define the sample, and 
no longer has a sharp cutoff. Figure |l] shows the original distribution of absolute magnitudes 
and the updated one. 

Column (11) — Radio (5 GHz)-to-optical {B band) flux ratio (Kellerman et al. 1989). 

Column (12) — Radio power at 5 GHz in W Hz~^ (Kellerman et al. 1989). 

Column (13) — H/5 full width at half maximum (FWHM), in km s^^ (Boroson & Green 
1992). 

Column ^i^;— ROSAT-PSPC X-ray spectral index a^ {F^ oc z/""-, Laor et al. 1997). 

Column ^ij;— [OIII] A5007 equivalent width, in A (EW, Boroson & Green 1992). 

Column (16)— Fell EW, in A (Boroson & Green 1992). 

Column (17) — [OIII] A5007 to H/? peak intensities ratio (Boroson & Green 1992). 

Column (18)— Rp EW, in A (Boroson & Green 1992). 

Column (19)—}iell A4686 EW, in A (Boroson & Green 1992). 

Column (20) — Optical (2500 A)-to- X-ray (2 keV) spectral index, Uox- We use our optical 
data to derive the 2500 A luminosity. X-ray data are taken from Tananbaum et al. (1986), 
Kriss (1988), and Laor et al. (1997). Laor et al.'s, based on ROSAT measurements, are the 
most reliable, and were preferred in case of disagreement. 

We have also studied the IR spectral indices (Neugebauer et al. 1987) and the polarization 
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levels (Berriinan et al. 1990) of the PG quasars. These parameters do not show significant 



correlations with any variability parameter we have measured (see § |3.3| ), so for brevity's sake 
they are not listed in Table 2. 

2.2 Observations and Data Reduction 

All objects were monitored for 7 years (1991 March through 1998 March) in the Johnson- 
Cousins B and R bands with the Wise Observatory 1 m telescope in Mitzpe-Ramon, Israel. 
The sampling properties (including the number of observing epochs, total period of obser- 
vations, and the median sampling interval) of each object are listed in Table 2. The typical 
(median over the entire sample) sampling interval is 39 days. 

Observations from 1991 March through 1994 January were obtained using an RCA 320 x 
520-pixel thinned CCD. Most of the data (1994 February through 1998 March) were obtained 
using a Tektronix 1024 x 1024-pixel back-illuminated CCD, with typical integration times of 
2-3 minutes in R and 3-4 minutes in B. During an intermediate period (1993 April through 
June), a TI 1000 x 1024-pixel CCD was used. Spectrophotometric monitoring of twenty-eight 
objects from our sample was carried out during the same period at Wise Observatory and at 
Steward Observatory (Maoz et al. 1994; Kaspi et al. 1996; Kaspi et al. 1999). Flux calibrated 
spectra for thirteen of the quasars are presently available. We multiplied those spectra with 
the Johnson-Cousins B and R transmission curves and integrated them over wavelength to 
derive B and R magnitudes, and thereby improved the temporal sampling of these objects. 

Owing to the large volume of data (~ 3500 CCD frames), reduction was carried out by 
specially written automatic routines that produced lightcurves from raw images with mini- 
mal interaction. Photometry was obtained with the point spread function (PSF) fitting rou- 
tines in the DAOPHOT (Stetson 1987) package within IRAF[j] applied to the bias-subtracted 
and fiat-field-corrected frames. This produced frame-dependent instrumental magnitudes for 
each quasars and for 5 — 10 reference stars in its field. One image of each quasar was used as 
a reference image. For each reference star, the difference between its magnitude in a given 
image and its magnitude in the reference image was added to the quasar magnitude in the 
image to obtain a quasar lightcurve relative to each reference star. These lightcurves where 

T IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, which are 
operated by AURA, Inc., under cooperative agreement with the National Science Foundation. 
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+65 50 38 


0.129 


64 


2273 


22 


14.9 


-23.8 


14.5 


-24.3 


-0.03 


23.3 


8450 


1.49 


20 


42 


1.99 


110 


2.5 


1.38 




16 17 56.2 


+ 17 31 34 


0.114 


56 


2308 


24 


15.5 


-22.8 


15.4 


-23.0 


-0.41 


22.8 


5330 


0.95 


7 


69 


0.48 


115 


5.1 


1.40 




16 26 51.5 


+55 29 05 


0.133 


30 


2308 


35 


15.7 


-23.1 


15.3 


-23.5 


-0.96 


22.1 


4490 


1.94 


9 


46 


0.56 


143 


11.8 


1.44 




17 00 13.4 


+51 53 37 


0.292 


54 


2403 


27 


15.4 


-25.2 


14.8 


-25.9 


0.38 


24.4 


2210 







78 


0.01 


55 


0.5 






17 04 03.5 


+60 48 31 


0.371 


40 


2447 


29 


15.6 


-25.5 


15.0 


-26.1 


2.79 


26.9 


6560 


0.96 


27 


5 


6.50 


28 


0.5 


1.45 




21 30 01.2 


+09 55 01 


0.061 


79 


2398 


16 


14.7 


-22.4 


14.3 


-22.8 


-0.49 


22.5 


2330 


1.95 


20 


70 


0.89 


109 


11.2 


1.29 




22 33 39.8 


+ 13 28 21 


0.325 


33 


2260 


42 


16.3 


-24.1 


16.4 


-24.0 


-0.55 


23.3 


1740 


1.73 


12 


60 


0.77 


67 


3.6 


1.49 




22 51 40.4 


+ 11 20 41 


0.323 


43 


2398 


29 


16.0 


-24.8 


15.5 


-25.3 


2.57 


26.4 


4160 




19 


26 


1.69 


82 


2.3 


1.66 





Table 2: Sample properties. See text for detailed description. 
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then averaged to obtain the final quasar hghtcurve. Reference stars whose magnitude differ- 
ence deviated by more than 3cr from the mean magnitude difference of the other reference 
stars at a certain epoch (presumably due to systematics such as bad pixels, cosmic rays, or 
variations of the reference stars) were excluded. Reference star measurements with photo- 
metric errors larger than 0.1 mag were also excluded. Photometric standards (Landolt 1992) 
were observed on several clear nights and were used to bring the lightcurves to an absolute 
scale. 

Lightcurves at B and R for each object are shown in Figure @. A machine-readable 
version of these data is available upon request from the authors. The uncertainty in each 
point is a combination in quadrature of the uncertainty of the PSF fitting error and the 
Icr scatter of the magnitude differences of the reference stars between that epoch and the 
reference epoch. Uncertainties are ~ 0.01 mag in the R band and somewhat higher (~ 0.02 
mag) in the B band, due to lower CCD sensitivity at the blue wavelengths. The uncertainties 
are dominated by the PSF fitting errors. Aperture photometry with large fixed apertures 
did not improve the accuracy. The uncertainty in the absolute photometric scale, which is 
typically 0.08 mag, is not included in the error bars. 

3 ANALYSIS 

3.1 Ensemble Variability Amplitude 

Figure ^ shows histograms of the brightness deviations about the mean levels in B, R, and 
B — R, along with their Icr scatter and the number of points contributing to the histogram. 
Also shown are the brightness deviation distributions of stars (one for each quasar field), Bs 
and Rs- The stellar distributions do not include spectrophotometric measurements (there is 
only one star per quasar in those measurements) and hence have fewer points. The standard 
deviations of the stellar distributions provide an independent estimate of the systematic 
error over all epochs for the entire sample: ass = 0.015 mag and a^s = 0.014 mag. Between 
February and November 1997 the telescope suffered from scattered-light problems as a result 
of a change in baffling. Data from this period have larger systematic errors (0.03 — 0.04 mag) 
due to the resulting flat-fielding problems. The following values can be regarded as the typical 
intrinsic variability amplitudes of the entire sample: 



^B = Jcrl - als = 0.14 mag (2) 



© 1998 RAS, MNRAS 000, 
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Figure 3. Distributions of the brightness deviations about the mean lightcurve levels of the B, R, and B — R measurements 
of the entire sample. Bg and Rg are the brightness deviation distributions of stars in the quasars' fields. 



^R = \^R - ^RS = 0-12 mag 



(3) 

ae-R = \j(yl-R - als - ^RS = 0-05 mag (4) 

In principle, the constant contribution from the host galaxies of the quasars may dilute their 
intrinsic variations. Bahcall et al. (1997) observed the host galaxies of ten quasars that are 
in our sample with the Hubble Space Telescope. From their measurements we estimate the 
galaxy contribution to be ~ 5% within our PSF fitting diameter of <, 3". The intrinsic quasar 
variations are therefore ~ 1.05 times larger than the numbers above. Nevertheless, the PG 
quasars in Bahcall et al. (1997) do not cover the entire range in redshift and luminosity 
of our sample, and so may not be completely representative. The host galaxies may also 
slightly change the colour of the quasars (see §P7^). 
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Another factor that can bias the variabihty parameters is the emission-hne contribution 
to the broad-band fluxes. In the PG sample, the typical emission-line contributions to the 
B and R fluxes are 5 — 10% and 10 — 15%, respectively. This includes both the constant 



narrow lines and the variable broad lines. We further discuss this issue in 5EI3 



3.2 Statistical Measures 

We define below the statistics we use to estimate the variability amplitude and time scale, 
and the method by which we test for the correlations between these parameters and other 
quasar parameters. 

Variability Amplitude - We have estimated the variability amplitude of each object 
in various ways. The standard deviation was calculated for each lightcurve, weighting each 
data point by its inverse variance. The maximal peak-to-peak difference between any two 
points in the lightcurve is another conventional estimator used by many authors (e.g., Netzer 
et al. 1996). Other meaningful measures are the variability mean and median, defined as the 
mean and median of all possible magnitude differences of a lightcurve. The median is a better 
estimator in terms of its lower sensitivity to outliers (Hook et al. 1994; Netzer et al. 1996). 
We have also calculated separately the medians, median{Am+) and median{Am^) , of the 
positive and the negative differences for each lightcurve, in order to test for asymmetries in 
the lightcurves. We define the net difference of the lightcurve by: 

Sm = median{Am^) — median{\Ara_\) . (5) 

The standard deviations and the mean, median, maximal, and net differences of the lightcurves 
are listed in Table 3. 

Variability Gradient - A simple estimate of variability is the variability gradient 
(Borgeest & Schramm 1994), defined as the magnitude change per unit time. Using all 
pairs of adjacent data points for each lightcurve (including the B — R lightcurves), we have 
calculated the weighted mean, median, and maximal gradient, in units of magnitude per 
year, in the quasars' rest-frame, and list them in Table 3. For the maximal gradient we 
took only pairs of points that are more than one day apart, since the combination of the 
small flux fluctuations (due to photometric uncertainties) over a few hours can translate into 
unrealistically large values in magnitudes per year. We note that the variability gradient is 
a crude estimator. Two completely different lightcurves in terms of typical time scale, such 
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as a periodic low-amplitude time series and a time series with a single high- amplitude flare, 
can give the same mean gradient. 

Power Spectrum Analysis - A time scale can also be measured, in principle, from 
features in a Fourier power spectrum. Previous observations of individual AGN, especially in 
the X-ray range, have yielded a featureless, power-law-shaped power spectrum, Pi^u) oc u^'^ , 
on time scales ranging from minutes to weeks (e.g.. Green et al. 1993; Paltani et al. 1997; 
Fiore et al. 1998). Recent study of NGC 3516 by Edelson & Nandra (1998) have shown that 
its power spectrum exhibits a break on time scale of ~ 30 days. We have performed a power 
spectrum analysis on the mean-subtracted lightcurves of the individual objects. 

There are algorithms (e.g., the Lomb-Scargle algorithm, Scargle 1989) designed for power 
spectrum analysis of unevenly-spaced data, whose basic approach is to place zeroes on an 
evenly-sampled grid at all points in the mean-subtracted lightcurve where there are no data. 
This approach works well when one is searching for periodicities. However, when dealing with 
noise-like spectra, such as those characterizing AGN, the above approach can distort the true 
slope of the power spectrum. This occurs because a lightcurve with zeroes where there are no 
data will be much "bluer" (i.e., with many short, large spikes, producing relatively too much 
power on short time scales) than it would be had it been evenly sampled. Stated differently, 
the long time scales do not get enough weighting in the Fourier sum, because the slow 
harmonics get multiplied mostly by zeroes. One can, alternatively, interpolate the lightcurve 
to an evenly-sampled grid, but this involves inventing artificially-smoothed variations in the 
observational gaps, with the effect of attenuating the high frequencies in the power spectrum. 
As a compromise between these two extremes (no interpolation and full interpolation) we 
have used a scheme (T. Mazeh 1998, private communication) where interpolated data are 
used only when the interval between the real and interpolated data points is smaller than a 
set fraction of the Fourier time scale being calculated. This appears to work well, in that only 
the desirable aspects of each approach are kept. Specifically, for every frequency v whose 
power is calculated, a linearly interpolated point is introduced in the lightcurve only if the 
time-difference. At, between the interpolated point and the closest real data point is smaller 
than l/2z/. Otherwise, the flux is set to zero at this epoch in the evenly-sampled grid. 

The two lowest frequency points of each power spectrum were ignored, due to insufficient 
statistics. The remaining points were binned logarithmically - the first bin was derived by 
averaging the two remaining lowest frequency points, the second by binning the next four 
points, eight points, etc. This binning scheme gives equal weightings to equal decades in 
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frequency (Edelson & Nandra 1998), when, e.g., fitting a power-law relation to the power 
spectrum. 

As an example, we show in Figure ^ three power spectra of PG0804+762. Figure ^ is 
the power spectrum resulting from the uninterpolated lightcurve. Due to the suppression of 
the low frequencies, it appears almost like a white-noise spectrum. Figure |^ is the power 
spectrum of the fully interpolated lightcurve. Here the higher frequencies are attenuated by 
the smoothing effect of the interpolation. Figure ^c shows the power spectrum resulting from 
our partial interpolation scheme. This power spectrum is similar to the "fully-interpolated" 
spectrum (Fig. ^) at low frequencies and to the "uninterpolated" spectrum (Fig. ||a) at 
high frequencies, with a smooth transition between them. In the lower panel we also show 
the power spectrum of the window function. It was obtained by sampling, with the same 
temporal pattern, white noise having our measurement errors and calculating the power 
spectrum with this algorithm. 

To test the reliability of our procedure in reproducing a given power-law index, we 
have carried out Monte-Carlo simulations. We have generated artificial lightcurves by sum- 
ming sine waves with amplitudes following a power-law relation and having random phases. 
Multiple lightcurves were produced by drawing different random phases. We then sampled 
each lightcurve with the actual temporal pattern of one of the observed quasars selected 
at random, and measured the distribution of the power-spectrum indices produced by our 
partial-interpolation procedure. We find from these simulations that the input index is best 
reproduced by fitting only frequencies < 0.01 day~^. The la scatter about the index that is 
input to the Monte-Carlo simulations is typically ±0.6. The mean of the Monte-Carlo dis- 
tribution of indices matches the input index when the denser sampling patterns (with more 
than 60 points) are used. For sparser sampling patterns the mean output index is systemat- 
ically different from the input index. The power spectrum of a daily sampled lightcurve in 
our simulations reproduced the input index with an accuracy of ~ ±0.1. 

We have applied this algorithm to our data in the frequency range of 0.001 day^^ < '^ < 
0.01 day~^. We find that the quasars generally have power spectra in this range of frequencies, 
that can be fitted with a single power-law with x^ values of unity order. Exceptions are 
PG1444±407, PG1415±451 {B band), and PG2233±134 {R band), whose power spectra 
show a break at ~ 300 days. The index values for the B and R lightcurves are in the ranges 
0.8 < 7_B < 3.3, and 0.6 < 7r < 3.1, and both have means of 7 = 2.0 ± 0.6. Based on 
the Monte-Carlo simulations, we therefore cannot distinguish the spectral indices of the 
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Figure 4. The R band power spectrum of PG0804+762, in logarithmic units: a) Power spectrum obtained without any 
interpolation of the unevenly sampled data, b) Power spectrum after full interpolation of the lightcurve to an evenly-sampled 
grid, c) Power spectrum resulting from the partial interpolation scheme described in the text. The power spectrum of the 
window function is shown for comparison in the lower panel. 
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individual PG quasars. However, we can conclude that the quasars have a limited range of 
a power-law indices, <, ±0.6. 

Autocorrelation Analysis - The autocorrelation function (ACF) is defined by 

ACF{t) = {{m{t) - (m)) ■ {m{t + r) - (m))), (6) 

where brackets denote a time average. The ACF measures the correlation of the lightcurve 
with itself, shifted in time, as a function of the time lag r. If there is an underlying signal 
in the lightcurve, with a typical time scale, then the width of the ACF peak near zero time- 
lag will be proportional to this time scale. Another function used in variability studies to 
estimate the variability at different time lags is the first-order structure function (SF; e.g., 
Trevese et al. 1994) defined by 

SF{r) ^ ^((m(t + r)-m(t))2). (7) 

There is a simple relation between the ACF and the SF, 

SF^{t) = 2{y - ACF{t)) , (8) 

where V is the variance of the lightcurve. We therefore perform only an ACF analysis on 
our lightcurves. 

A "correlation time scale" can be defined by the width of the ACF peak. For a power- 
law-type power spectrum, that does not have a characteristic time scale, the power-law index 
and the S/N determine the correlation time scale. In our analysis we choose the zero-crossing 
time of the ACF as the correlation time scale, which designates the time when the signal 
has entirely "forgotten" its past. We have estimated the autocorrelation function (ACF) 
using the z-transformed discrete correlation function (ZDCF, Alexander 1997). Alexander 
(1997) has shown that this method is statistically robust even when applied to very sparsely 
and irregularly sampled lightcurves {nobs > 12). The ZDCF was calculated for all of the 
lightcurves. We then used a least-squares procedure to fit a fifth-order polynomial to the 
ZDCF, with the constraint that ACF{t = 0) = 1. We used the ZDCF fit to evaluate the 
zero-crossing time in the observer's frame. As an example of the autocorrelation analysis, 
the ZDCF with its polynomial fit for PG0804-I-762 is shown in Figure ^. ACF time scales 
are listed in Table 3. 



The statistics described above are listed in Table 3 in the following format: 
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Figure 5. The ZDCF (discrete points) and its fit (fifth order polynomial least-squares fit, solid line) for the PG0804+762 R 
band data. 

Column (1) — PG name (shortened). 

Column (2) — Standard deviation of the hghtcurves in the B band. 

Columns (3)- (5) — Mean, median, and maximal differences of the hghtcurves in the B 
band. 

Columns (6)- (8) — Mean, median, and maximal rest-frame gradients of the hghtcurves 
in the B band. 

Columns (9)- (15) — Same as columns (2)- (8), in the R band. 

Columns (16) -(22) — Same as columns (2)- (8), for the B — R colour. 

Column (23) — Net difference of the hghtcurves, in the B band. 

Column (24) — Same as column (23), in the R band. 

Column (25) — Observer's frame ACF time scale, in the B band, in days. A 'V is 
marked for objects whose ACF did not reach zero correlation. These quasars were excluded 
from the correlation analysis. 

Column (26) — Same as column (25), in the R band. 

Correlation Analysis - We have tested for the existence of correlations among all of 
the measured quantities described above. We have also searched for correlations with other 
parameters known for the PG sample (Kellerman et al. 1989; Boroson & Green 1992; Laor et 
al. 1997) that are listed in Table 2. Two additional data sets were included in the correlation 
analysis - the IR spectral indices from Neugebauer et al. (1987), and the polarization data 
from Berriman et al. (1990). Of special interest are the X-ray spectral index {ax) and the H/3 
FWHM, that were found by Fiore et al. (1998) to be correlated with the X-ray variability 
amplitude on short (2-20 days) time scales, but uncorrelated on longer time scales. We also 
sought correlations in the subclasses of the radio-loud quasars (RLQ, quasars with radio-to- 
optical flux ratio, R > 10) and the radio-quiet quasars (RQQ, R < 10) in our sample. There 
are 7 RLQ and 35 RQQ in our sample. 
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0.07 
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Table 3: Variability properties of the lightcurves. See text for more details. 
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We chose Spearman's rank-correlation coefficient over Pearson's correlation coefficient 
because it tests for a general monotonic relation rather than only a linear one. However, in 
cases of bimodality Spearman's coefficient may fail to detect correlations (e.g., the variability- 
amplitude/luminosity correlation which we discuss below). Table 4 shows the probability 
matrix of the quasar properties, where each entry gives the probability that the data are 
drawn from two uncorrelated populations. 

There are 322 correlations in Table 4, so naively one expects 16 spurious correlations with 
Pr < 5%. We have found 51 correlations with probabilities < 5%. However, only 11 of the 
correlations are between independent parameters. The remaining 40 are between dependent 
parameters {apriori known to be correlated) or degenerate parameters (functions of each 
other). There are only two independent variability parameters - the variability amplitude and 
the variability time scale. This is because the variability amplitude parameters for each band 
are degenerate variables, since they are functions (although complicated) of each other, and 
the two bands are dependent. Therefore, all of the variability amplitude parameters should 
be counted only once. The variability time scales in the two bands are also dependent. 
Similarly, the luminosities in both bands are dependent variables, and both are coupled to 
the redshift, since the PG sample is fiux-limited (see below); the radio-loudness and the radio 
power are both coupled to the luminosity (Hooper et al. 1996, see below); the H/5 FWHM 
and a^ are also correlated (Laor et al. 1997), and ax and Oox are dependent variables, too. 
Therefore, there are 8 independent spectral properties. This gives a total of 16 independent 
correlations, for which 1.6, and 0.8 are expected to be spurious for Pj. < 10% and Pr < 5%, 
respectively. Thus, given that we find 11 independent correlations, we expect < 1.1 of them 



to be spurious. We discuss the significant correlations listed in Table 4 in § p73 . 

We have also tested for correlations between the change in the B — R colour and the 
change in brightness in individual objects. It is established in the best-studied Seyfert galax- 
ies that as the continuum source brightens, its UV-optical spectrum becomes harder, i.e., the 
variability amplitude increases with the emitted frequency (e.g., Maoz et al. 1993, Romano 
& Peterson 1998). However, we have found only one work, which had demonstrated this 
behaviour in a sample of quasars (Cutri et al. 1985). We have calculated the correlations 
between the B — R colour and both the B and the R magnitudes for each object. The results 
are shown in Table ^. 
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Table 4: Spearman's rank order correlation coefficient matrix for various quasar properties. Shown are only the probabilities for random 
correlations, in percents. Probabilities lower than 5% are highlighted. A minus sign represents an anticorrelation. The columns are the 
spectral properties and the rows are the variability properties, and are marked as in Tables 2 and 3. Lb and Lr are the luminosities in the 
rest-frame B and R bands, respectively. See text for a detailed description. 
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Table 5. Spearman's correlation coefficients and probabilities for tiic B — R colour changes vs. the B and R magnitude changes 
for individual objects. 



3.3 Correlation Results 

In this section we review the most significant correlations detected in the data, as listed in 
Tables 4 and ^ The section is divided on the basis of the main variability properties measured 
for the sample - the correlation of changes in colour with changes in magnitude, variability 
amplitude correlations, variability time scale correlations, variability gradient correlations, 
and colour variability correlations. All of the probabilities referred to in this section are 
taken from these tables, unless otherwise stated. 

Colour-Change Correlations - In Table |^ we have listed Spearman's correlation coef- 
ficients and probabilities for the correlations between the colour changes and the magnitude 
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Figure 6. Deviations from the lightcurve means of the B — R colour vs. the deviations of the B and R magnitudes, for the 
data points of the 21 quasars that become bluer as they brighten. The number of data points and the best fit slopes for each 
band are specified. 

changes in the B and R bands for individual objects. Thirty eight objects in Table ^ show 
significant correlations {Pr < 5%) of the colour changes with the changes in the B mag- 
nitude, while only twenty three objects display this behaviour in the R band. We suspect 
that the apparently-significant A{B — R) vs. AB correlations of some of the quasars are due 
to the correlated errors between A{B — R) and AB, while some significant A{B — R) vs. 
AR correlations are likely lost due to the corresponding anticorrelated errors (note that PG 
1229-^204, PG 1444+407, and PG 1700+518 become redder when they brighten, probably 
due to this effect, because it is seen only in the R band). In any case, it is clear that a 
large fraction of the PG quasars become bluer as they brighten. In order to quantify this 
behaviour, we have divided the sample into a subgroup (1) of the 21 quasars which get bluer 
as they brighten (-Pr(-R) < 5.5%), and a subgroup (2) of the 21 quasars which do not. Figure 
P shows the deviations from the lightcurve means of the B — R colour vs. the deviations of 
the B and the R magnitudes, for the data points of subgroup 1. This figure manifests very 
significant correlations. The best linear fits are marked in the figure. 

To investigate if any of the quasar parameters is driving this division of our sample, we 
checked for correlations between the various quasar parameters, referred to in the preceding 
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Figure 7. Distributions of the median amplitude difference in the B band. Empty bars are objects which become bluer when 
they brighten, and filled bars are objects which do not exhibit a correlation between their changes in colour and brightness. 

sections, and the probabilities Pr{R) in Table ^. We find that the quasars with the larger 
variability amplitudes are those that tend to become bluer when brightening in R. Figure |^ 
shows the histograms of the median difference in the B band, which represents the variability 
amplitude, with the two subgroups distinguished. A Student's t-test on the means of the two 
subgroups shows them to be different at a significance level > 99%. From the histograms it 
is evident that for AB ^ 0.2 mag, all objects become significantly bluer when they brighten. 
The range 0.1 ^ AB <^ 0.2 mag is occupied by both populations. For AB ^0.1 mag, the 
uncertainty of the colour measurements becomes dominant [A{B — R) = (0.266 ±0. 007) Ai?, 
Figure |^, and a definitive statement about the colour variations cannot be made. Some 
role in this effect may be played by the colours and magnitudes of the host galaxies, but 
investigating this will require improved imaging data for the hosts of this sample. 

Variability Amplitude Correlations - We find that the variability amplitude in its 
various forms is weakly anticorrelated with the source optical luminosity with probabilities 
of Pr = 9 — 20% for the R band and of P,. = 13 — 26% for the B band. This is in accord with 
previous claims by most of the authors mentioned in Table [^ (Pica & Smith 1983; Trevese 
et al. 1994; Hook et al. 1994; Hawkins 1996; Cristiani et al. 1997). The significance of this 
correlation is weaker than that found in the study of Pica & Smith (1983) - Pj. = 7% for 
RQQ, and weaker than found by Hook et al. (1994) and Cristiani et al. (1997), which is 
Pr < 1%. The same weak trend exists also between the variability amplitude and the radio 
power (at 5 GHz), at a significance level of P,- = 10 — 11% (However, Impey et al. (1991) have 
found that for RQQ the variability amplitude is correlated with the radio power). In Figure 



© 1998 RAS, MNRAS 000, 



30 Giveon et al. 



0.1 



0.0 






o 
o o 



0.1 



0.0 



••• 

• •• 



o 
o 
o 



29.25 29.75 30.25^ 3^.75 31.25 20 

logL„ [erg s" Hz" ] 



22 24 26 

logPj [W Hz ] 



0.2 



0.0 











• 










• • 








U.2 


• 

•• • 

• • 








n.1 


•ttJ; 






Q3o 




•• 




o 




• • • 




O 






1.0 



1.2 
1+Z 



1.3 



1 

logRj 



Figure 8. R band standard deviation vs. rest-frame R band luminosity (Ln), redshift (1 + z), total radio power at 5 GHz 
(Ps), and the radio-Ioudness (Rs). Filled circles are radio-quiet objects and empty circles are radio-loud objects. 

H we show the R band standard deviation vs. luminosity, redshift, radio power, and radio- 
loudness. The lower mean variability amplitude of the RLQ is clear, even though there is no 
significant monotonic trend according to Spearman's rank test {Pr = 14—24%). However, the 
quasars with the highest optical luminosity in the PG sample are RLQ (as already noted by 
Hooper et al. 1996; our data give a correlation of r^ = 0.512, P^ = 5.3 x 10~^ between Lji and 
R), so this may be a luminosity effect. From the a^ vs. logLn plot, no significant difference 
exists in variability amphtude between objects in the RLQ and the RQQ subgroups that have 
comparable luminosities (ctrl = 0.08 ±0.03 and ctrq = 0.10 ±0.04, and are indistinguishable 
at the 80% significance level according to Student's t-test). Hence we cannot say whether the 
low amplitude of the RLQ is the result of their high optical luminosity or their radio power, 
and there is no evidence for the RLQ having larger variability amplitude than the RQQ, in 
agreement with the findings of Pica & Smith (1983) and Borgeest & Schramm (1994). 

We find no significant correlation between the variability amplitude parameters and the 
source redshift. However (see Figure |]), here too there may be a bimodality, with z < 0.2 
objects having larger aR than z > 0.2 objects {<Jz<o.2 = 0.12 ±0.05 and crz>o.2 = 0.09 ±0.04, 
which are different at the 95% significance level according to Student's t-test). We note 
also that in the lower redshift bin {z < 0.2) the trend of the variability amplitude-redshift 
correlation is opposite to the trend in the higher redshift bin {z > 0.2). The PG sample 
suffers from a well-known luminosity-redshift bias, as do other flux-limited samples (Schmidt 
& Green 1983). Our data give a correlation coefficient r^ = 0.608, Pr = 1.9 x 10"^, between 
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Figure 9. B band net difference histogram. Dashed-lined bars correspond to the four objects (PG 0026+129, PG 1048+342, 
PG 1427+480, and PG 1626+554) exhibiting long-term monotonic brightening. 

the rest-frame R band luminosity and the redshift. Thus, the correlations of variability 
amplitude with luminosity and redshift cannot be disentangled in our sample. 

The variability amplitude is anticorrelated with the net difference (equation ^) with 
significance ^ 99% (not shown in Table 4). This means that objects that show larger varia- 
tions exhibit asymmetry (negative 5B) in their variations - on average the brightening takes 
longer then the fading. However, this trend is dominated by four objects (PG 0026+129, PG 
1048+342, PG 1427+480, and PG 1626+554) whose lightcurves exhibit long term mono- 
tonic variations and therefore have large 5B and 5R. Excluding them from the correlation 
analysis reduces the above significance to 95 — 99%. The net difference histogram is shown 
in Figure ^, with the four monotonically rising objects highlighted. When these four objects 
are excluded, the mean of the distribution is 5B = —0.02, which is consistent with zero 
according to Student's t-test. 

A strong result of our study is a correlation of the variability amplitude with the H/3 EW. 
It emerges from all of our variability amplitude test parameters, and in both bands (see Table 



4). The average difference in the R band is plotted vs. H/5 EW in Figure |T0|. A possible, but 
unlikely, explanation for this trend is the emission-line contribution to the broad-band fluxes. 
In objects with large emission-line EW, small intrinsic optical continuum variability may be 
accompanied by larger variations in the ionizing continuum, producing large variations in the 
optical emission lines (after a delay of several months, due to the light-travel time across the 
broad-line region; Kaspi et al. 1996). The delayed and smoothed response could also explain 
the larger ACF time scales in quasars with larger EW (see below) . However, in all AGN that 
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Figure 10. The median difference in the R band vs. H/3 EW. Filled circles are radio-quiet objects and empty circles are 
radio-loud objects. 

have been monitored spectroscopically to date, the fractional amphtude of the emission-hne 
variations is always considerably smaller than that of the optical continuum variations, once 
known constant components such as narrow lines and galaxy light have been accounted for 
(e.g., Maoz 1997). Spectrophotometry of several of the PG quasars suggests this behaviour 
applies to this sample as well (Maoz et al. 1993; Kaspi et al. 1996). Spectrophotometric 
monitoring data for a subsample of the quasars may allow us to understand this correlation 
(Kaspi et al. 1999). The correlation with the variability amplitude is also detected for Hell 
A4686 EW, [OIII] A5007 EW, and [OIII] A5007-to-H/3 peak intensities ratio {B band only), 
though at weaker significance levels (see Table 4). Apparently, this trend is not a combination 
of the variability amplitude-luminosity anticorrelation and the Baldwin effect (EW oc L~^ 
for various lines; Baldwin 1977, Osmer et al. 1988), since H/? is one line that does not exhibit 
the Baldwin effect (Binette et al. 1993). 

Variability Time Scale Correlations - We find that the rest-frame ACF variability 
time scale is correlated with the source luminosity, but not correlated with redshift. The 
correlation with the luminosity is opposite to the ones found by Trevese et al. (1994) and 
Netzer et al. (1996). The strongest correlation (P^ = 2.1%) is between the rest-frame ACF 
time scale in the R band and the rest-frame B band luminosity, which are plotted in Figure 
n]. Four quasars (all RQQ), whose ACF did not reach zero correlation (see Table 3), were 
excluded from the calculations. 

A weak anticorrelation exists between the B band rest-frame ACF time scale and the 
X-ray spectral index, a^ {Pr = 7.8%), shown in Figure |12. 



Variability Gradient Correlations - The median and the mean rest-frame gradients 



are anticorrelated with luminosity. This is shown graphically in Figure |T3. 
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Figure 11. Rest-frame ACF time scale in the R band vs. the rest-frame B band luminosity. Filled symbols are radio-quiet 
objects and empty symbols are radio-loud objects. 
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Figure 12. Rest-frame ACF time scale in the B band vs. o^. Filled symbols are radio-quiet objects and empty symbols are 
radio-loud objects. 

We detect a few significant correlations between the rest-frame variability gradient and 
a few parameters, which seem peculiar, in that related measures, e.g., the mean and median 
gradients, show no comparable correlations. For example, the mean R band gradient and 
the [OIII] A5007 to H/3 peak intensities ratio are significantly correlated (P^ = 4.0%) but 
the median R band gradient is not correlated with that parameter. 

Colour Variability Correlations - We have also tested for correlations of the quasar 
parameters with the B — R colour variation properties (see Tables 3 and 4). All of the 
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Figure 13. Mean rest-frame gradient in magnitudes per year vs. monochromatic luminosity, Lb.. Filled circles are radio-quiet 
objects and empty circles are radio-loud objects. 
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Figure 14. Relative mis amplitude of the B — R colour vs. rest frame R band luminosity, L^. Filled circles are radio-quiet 
objects and empty circles are radio-loud objects. 

parameters used to estimate the colour variability amplitude show significant anticorrela- 
tions with the optical luminosities, and some of them show anticorrelations with the radio 
power and the radio-loudness. The anticorrelations with the radio power may be due to the 
radio-loudness-luminosity trend combined with the colour-luminosity anticorrelation which 
is exhibited by the sample (Mushotzky & Wandel 1989). In Figure |1^ we show the B — R 
rms amplitude vs. the R band rest-frame luminosity [Pj. = 1.3%). 

Less luminous objects tend to have larger colour gradients, as is evident from the colour 
mean gradient-luminosity anticorrelation (P^ < 2%) seen in Table 4. The maximal colour 
gradient is correlated with the radio-loudness {Pr = 3.0%) but not with the optical luminos- 
ity. These correlations suggest that objects that are more radio-loud have colour lightcurves 
with flickering activity, but the total variance of their colour lightcurves is lower. 

We have checked for correlations between the variability parameters and two additional 
data-sets - the IR spectral indices (Neugebauer et al. 1987) and the polarization levels 
(Berriman et al. 1990) of the PG quasars, but no significant correlations were found. The 
PG quasars have low polarizations and very low polarization variability. The polarization 
is mostly due to scattering, as opposed to BL Lac objects, where synchrotron emission is 
responsible for it. 

3.4 Comparison With Other Studies 

As already outlined in §|I], previous studies of quasar variability have often reached conflicting 
results. For example, among the fourteen studies listed in Table |T|, six found an anticorre- 
lation between the variability amplitude and the source luminosity, and five did not detect 
such a correlation. The relation between the variability amplitude and the source redshift is 
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even more confounding - there are three positive correlations, one anticorrelation, and seven 
studies that did not detect any correlation. 

In the present work we have detected a weak anticorrelation between the variability 
amplitude and the source luminosity, which agrees with most of the studies listed in Table |1|. 
Taking into account only studies of samples with luminosities and redshifts comparable to 
ours (Bonoli et al. 1979; Giallongo et al. 1991; Cimatti et al. 1993; Trevese et al. 1994), 
our work agrees with the minority of Trevese et al. (1994). No significant correlation was 
detected in our data between the variability amplitude and the source redshift, in accord 
with most of the previous studies. However, the sources redshifts of the our sample are in 
a rather narrow range {z < 0.4) compared with other studies. In any case, we find that 
these trends are weak and have intrinsically large scatter, as evidenced by the fact that they 
remain weak even when the observational noise is minimized. 

4 CONCLUSIONS 

We have presented the results of a seven-year monitoring program of a sub-sample of 42 
quasars from the PG sample, and studied their variability properties in the B and R bands. 
Our main findings are: 

1. All 42 objects have varied during the seven years with 5% < ctb < 34% and 4% < 
an < 26%. The rms amplitudes and colour variation for the entire sample are as = 14%, 
aji = 12%, and ctb-r = 5%. 

2. The power spectra of most of the objects have a power-law shape {Pi, oc u^"^) on time 
scales of 100 to 1000 days. The mean and rms values for the entire sample are 7 = 2.0 ± 0.6, 
and this spread is comparable to our uncertainties in determining 7. 

3. Our data reproduce the previously-reported anticorrelation between the variability am- 
plitude and the luminosity, but this trend is weak. A weak anticorrelation between the vari- 
ability amplitude and the redshift is also detected but may be the result of the luminosity- 
redshift bias of the PG sample. 

4. Larger variations tend to be mildly asymmetric - the brightening phase is longer than 
the fading phase. 

5. The ACF time scale is correlated with luminosity. More luminous quasars thus exhibit 
variations with longer time scales. 

6. The spectra of about half of the quasars become harder (bluer) when they brighten, 
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and the colour changes of the rest are not correlated with brightness changes. The quasars 
that become bluer tend to be those that are most variable. 

7. The RLQ tend to have lower variability amplitudes than the RQQ, but this may be due 
to the fact that the RLQ in this sample are more luminous than the RQQ. 

8. Quasars that are more luminous and more radio- loud have smaller colour variations on 
average, but RLQ tend to exhibit flickering in their B — R lightcurves. 

9. The median rest-frame variability gradient is anticorrelated with luminosity and with 
redshift. In other words, variations are less abrupt at higher luminosities and higher redshift. 
We are not able to disentangle these anticorrelations from the luminosity-redshift bias. 

10. Objects with larger equivalent widths of the H/3, [OIII] A5007, Hell A4686 emission lines 
have larger variability amplitudes. 

The overall picture emerging from these results is that the spectral and variability pa- 
rameters which describe the PG quasars have two extremes. The higher-luminosity quasars, 
which are also radio- loud, bluer, and at higher redshift, exhibit lower amplitudes of varia- 
tions with longer time scales, and their colour lightcurves exhibit flickering activity. At the 
other extreme, the lower-luminosity quasars, which in the PG sample are radio-quiet and at 
lower redshift, have a wide range of variability amphtudes and time scales but, on average, 
vary with larger amplitudes. 

We point out, however, that the trends we have identified are weak and the scatter is 
always much greater than the trend itself. This holds even for the strongest trend we have 
found, between variation amplitude and emission-line EW. The only trend with relatively 
little scatter is the tendency of about half of the quasars to become bluer when they brighten 
and redder when they fade. 

A main conclusion of this work is that none of the quasar parameters that we have 
examined are very reliable predictors of variability properties, even when observational noise 
(photometric errors, sampling frequency, and sample selection) are low. Future studies of 
other samples and other parameters may reveal improved variability predictors. 

The mechanism behind quasar variability is unknown. Some models have been proposed 
to explain it, but none of them comprehensively account for the variability phenomenon and 
the entire range of other quasars properties. Currently debated models include models, which 
suggest that the variability is caused by instabilities in an accretion disk (e.g., Mineshige 
et al. 1994a,b; Bao & Abramowicz 1996; Kawaguchi et al. 1998), starburst models, which 
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argue that quasar variability is formed by superposing the independent signals from many 
supernova explosions (e.g., Cid Fernandes et al. 1996; Aretxaga et al. 1997), and microlensing 
models, which posit that much of the observed variability could be caused by external agents, 
such as intervening lenses (e.g., Hawkins 1996; see however, Alexander 1995). We defer to 
a future paper a detailed comparison of our results to the predictions of these and other 
models. 

ACKNOWLEDGEMENTS 

We thank Buell Jannuzi for contributing to the observations at Steward Observatory, and 
for valuable comments. Tal Alexander, Ari Laor, Tsevi Mazeh, Rudy Schild, T. Wang, and 
the anonymous referee are acknowledged for their useful input. We are grateful to Sammy 
Ben-Guigui, John Dan, Peter Ibbetson, Ezra Mashal, and Haim Mendelson of the Wise 
Observatory staff, and to Anna Heller for their assistance with the observations. This work 
was supported by grants from the Israel Science Foundation. 

REFERENCES 

Alexander, T. 1995, 274, 909 

Alexander, T. 1997, in "Astronomical Time Series" Eds: D., Maoz, A., Sternberg, & E.M., Leibowich (Dordrecht: Kluwer), 163 

Allen, C.W. 1983, "Astrophysical Quantities" (London: Athlone) 

Aretxaga, I., Cid Fernandes Jr., R., Terlevich, R.J. 1997, MNRAS, 286, 271 

Bao, G. & Abramowicz, M.A. 1996, ApJ, 465, 646 

Bahcall, J.N., Kirhakos, S., Saxe, D.H., & Schneider, D.P. 1997, ApJ, 479, 642 

Baldwin, J.A. 1977, ApJ, 214, 679 

Berriman, G., Schmidt, G.D., West, S.C, & Stockman, H.S. 1990, ApJS, 74, 869 

Bershady, M.A., Trevese, D., & Kron, R.G. 1998, ApJ, 496, 103 

Binette, L., Fosbury, R.A., & Parker, D. 1993, PASP, 105, 1150 

Bonoli, F, Braccesi, A., Federici, L., Zitelli, V., & Formiggini, L. 1979, A&AS , 35, 391 

Borgeest, U. & Schramm, K.J. 1994, A&A, 284, 764 

Boroson, T.A. & Green, R.F. 1992, ApJS, 80, 109 

Bregman, J. 1990 A&A Rev., 2, 125 

Cid Fernandes Jr., R., Aretxaga, I., Terlevich, R.J. 1996, MNRAS, 282, 1191 

Cimatti, A., Zamorani, G., & Marano, B. 1993, MNRAS, 263, 236 

Clavel, J. et al. 1991, ApJ, 366, 64 

Cristiani, S., Trcntini, S., La Franca, L., & Andreani, P. 1997, A&A, 321, 123 

Cutri, R.M., Wisniewski, W.Z., Rieke, G.H., & Lebofsky, M.J. 1985, ApJ, 296, 423 

Edelson, R. & Nandra, K. 1998, preprint 

Fiore, F., Laor, A., Elvis, M., Nicastro, F., & Giallongo, E. 1998, ApJ, 503, 607 

Giallongo, E., Trevese, D., & Vagnetti, F. 1991, ApJ, 377, 345 

© 1998 RAS, MNRAS 000, 



38 Giveon et al. 

Goldschmidt, P., Miller L., La Franca F., & Cristiani, S. 1992, MNRAS, 256, 65 

Green, A.R., McHardy, I.M. & Lehto, H.J. 1993 MNRAS, 265, 664 

Hawkins, M.R.S. 1996, MNRAS, 278, 787 

Hook, I.M., McMahon, R.G., Boyle, B.J., & Irwin, M.J. 1994, MNRAS, 268, 305 

Hooper, E.J., Impey, CD., Foltz, C.B., & Hewett, P.C. 1996, ApJ, 473, 746 

Impey, CD., Lawrence, C.R., & Tapia, S. 1991, ApJ, 375, 46 

Kaspi, S., Smith, P.S., Maoz, D., Netzer, H., & Jannuzi, B.T. 1996, ApJ, 471, L75 

Kaspi, S. et al. 1999, in preparation 

Kawaguchi, T., Mineshige, S., Umemure, M., & Turner, E.L. 1998, ApJ, 504, 671 

Kellermann, K.I., Sramek, R., Schmidt, M., Shaffer, D.B. & Green, R. 1989, ApJ, 98, 1195 

Kohler, T., Groote, D., Reimers, D., & Wisotzki, L. 1997, A&A, 325, 502 

Kriss, G.A. 1988, ApJ, 324, 809 

Landolt, A.U. 1992, AJ, 104, 340 

Laor, A., Fiore, F., Elvis, M., Wilkes, B., & McDowell, J.C 1997, ApJ, 477, 93 

Maoz, D. et al. 1993, ApJ, 404, 576 

Maoz, D., Smith, P.S., Jannuzi, B.T., Kaspi, S., & Netzer, H. 1994, ApJ, 421, 34 

Maoz, D. 1997, in "Emission Lines in Active Galaxies" Eds: B.M. Peterson, F.-Z., Cheng, & A.S., Wilson (PASP), p. 138 

Miller, P., Rawlings, S., Saunders, R., & Eales, S. 1992, MNRAS, 254, 93 

Mineshige, S., Ouchi, N., & Nishimoti, H. 1994a, PAS J, 46, 97 

Mineshige, S., Takeuchi, M., & Nishimoti, H. 1994b, ApJL, 435, 125 

Mushotzky R.F. & Wandel, A. 1989, ApJ, 339, 674 

Netzer, H. & Sheffer, Y. 1983, MNRAS, 203, 935 

Netzer, H. et al. 1996, MNRAS, 279, 429 

Neugebauer, G., Green, R.F., Matthews, K., Schmidt, M., Soifer, B.T., & Bennet, J. 1987, ApJS, 63, 615 

Osmer, P.S., Porter, A.C, Green R.F. 1988, BAAS , 20, 968 

Paltani, S., Courvoisier, T.J.L., Blecha, A., & Bratschi, P. 1997, A&A, 327, 539 

Pica, A.J. & Smith, A.G. 1983, ApJ, 272, 11 



Romano, P. & Peterson, B.M. 1998 (submitted, |astro-ph/980619C| ) 

Scargle, J.D. 1989, ApJ, 343, 874 

Schmidt, M. & Green, R.E. 1983, ApJ, 269, 352 

Siemiginowska, A. & Elvis, M. 1997, ApJ, 482, L9 

Steston, P.B. 1987, PASP, 99, 191 

Tananbaum, H., Avni, Y., Green, R.F., Schmidt, M., & Zamorani, G. 1986, 305, 57 

Trevese, D., Kron, R.G., Majewski, S.R., Bershady M.A., & Koo, D.C 1994, ApJ, 433, 494 

Uomoto, A.K., Wills, B.J., & Wills, D. 1976, AJ, 81, 905 

Usher, P.D. 1978, ApJ, 222, 40 

This paper has been produced using the Royal Astronomical Society/Blackwell Science M^iK 
style file. 



© 1998 RAS, MNRAS 000, 



